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We derive the effective action describing the long- wavelength low-energy collective modes of quasi- 
one-dimensional spin-density-wave (SDW) systems, starting from the Hubbard model within weak 
coupling approximation. The effective action for the spin-wave mode corresponds to an anisotropic 
non-linear sigma model together with a Berry phase term. We compute the spin stiffness and the 
spin-wave velocity. We also obtain the effective action for the sliding mode (phason) taking into 
account the density fluctuations from the outset and in presence of a weak external electromagnetic 
field. This leads to coupled equations for the phase of the SDW condensate and the charge density 
fluctuations. We also calculate the conductivity and the density-density correlation function. 
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I. INTRODUCTION 



Materials with highly anisotropic crystal structure 
often exhibit, density-wave (DW) instabilities at low 
temperature Jll While inorganic linear chain compounds 
usually develop a charge-density-wave (CDW) instabil- 
ity, several group of organic conductors present a SDW 
ground state. Well known examples include transition 
metal bronzes such as K0.3M0O3 for CDW systems, 
and the Bechgaard salts (TMTSF)2X (X=C104,PF6) for 
SDW systems. 

In quasi-one-dimensional (quasi-lD) systems, DW in- 
stabilities result from nesting properties of the Fermi 
surface. Indeed, in this case the particle-hole response 
function exhibits a logarithmic singularity ~ \ti{Eq/T), 
where Eq is an ultra-violet cutoff of the order of the 
bandwidth. In presence of repulsive electron-electron or 
electron-plionon interactions, this leads to an instability 
of the metallic phase at low temperature. CDW ground 
states resulting from elcctron-phonon mteraction were 
first discussed by Frohlichnl and Pcicrls,0 while the pos- 
sibility of a SDW ground state due to repulsive electron-, 
electron interaction was first postulated by Overhauser.Q 

In the SDW ground state, quasi-particle excitations 
exhibit a gap 2Ao, where Aq is the SDW order parame- 
ter. The low-energy behavior of the system is then dom- 
inated by low-lying collective modes. In presence of an 
incommensurate SDW, two continuous symmetries are 
spontaneously broken: the translational symmetry and 
the rotational symmetry in spin space. This leads to the 
existence of two gapless Goldstone modes: a sliding mode 
(phason) and a spin- wave mode (magnon). Contrary to 
the case of superconductors, collective modes in DW sys- 
tems directly couple to external fields, so that they can 
easily be observed in various experiments. For instance, 
the non-linear dc conductivity is a manifestation of the 
existence of a phason mode which is pinned by impurities 
in real systems. 



In CDW systems, collectiiie modes were first studied 
by Lee, Rice and Andersonjj From the computation of 
the Green functions, they deduced the existence of a 
gapped amplitude mode and a gapless phase mode (pha- 
son) . The latter corresponds to a sliding of the CDW and 
would lead to an infinite conductivity in a clean systerp. 
(i.e. with no impurity) as first proposed by Frohlich.H 
Fukuyama then proposed an effective phase Lagrangian 
C{6) determining the dynamics of the phason J3 where the 
condensate phase 9 determines the position of the DW 
with respect to the underlying crystal lattice. The first 
attempt to derive CiO) rigorously is due to Brazovskii 
and Dzyaloshinskii.Ll The phase Lagrangian was latter 
used to study the interaction of the CDW with impuri- 
ties and the mechanism of pinning/depinning which is at 
the origin of the non-linear dc conductivity observed in 
transport experimentsJj'EI 

To a large extent, the analysis of collective mc 
CDW systems can be transposed to SDW systemsBEj" 
The amplitude mode has a gap 2Ao and is therefore 
strongly damped due to the coupling jjjith quasi-particle 
excitations above the mean-field gap.EJ As in CDW sys- 
tems, the phason corresponds to a sliding of the DW 
and leads to an infinite Frohlich conductivity in the ab- 
sence of impurities. However, since the SDW instability 
is driven by electron-electron interactions, the coupling to 
the lattice plays no role and its mass is not renormalized 
by phonons. The pinning by impurities is also weaker, 
since the SDW couples to charge inhomogeneities only 
to second order. The spin-wave mode, which is specific 
to SDW system_iS|-also obtained from the poles of the 
Green functions.EJdiS 

The functional integral formalism has also gmveii use- 
ful to study collective modes in DW systems.LliiZrtj The 
main advantage of this formalism is that it allows to de- 
rive the effective Lagrangian of the Goldstone modes from 
first principles. Both the phase Lagrangian C{6) and the 
effective Lagrangian for the spin-wave mode can be ob- 



tained within this formahsm. 

In this paper we derive the effective action describing 
the low-energy collective modes of quasi-lD SDW sys- 
tems, starting from the Hubbard model with weak on-site 
interaction. This effective action describes the behavior 
of the system at energies much smaller than the mean- 
field gap Ao or, equivalently, at length scales much larger 
than the mean-field coherence lengths vp / ^o and wj_/Ao 
{vp and v±^ are the velocities along and across the con- 
ducting chains). The improvement with respect to pre- 
vious works is twofold. First, we show that the effective 
action for the spin-wave mode is given by an anisotropic 
non-linear sigma model (NLuM) together with a Berry 
phase term. Suci-a result is known for the isotropic 2D 
Hubbard model,Ej but has not been derived for weakly 
coupled chain systems. Instead, it is generally assumed 
that the spin dynamics cjan be deduced from an effective 
Heisenberg Hamiltonian.EJ Second, we introduce the long- 
wavelength charge-density field from the outset. We thus 
obtain an effective Lagrangian C(9,p) which is a func- 
tional of two independent fields: 0, which is the phase 
of the SDW condensate and the charge-density field p. 
This yields coupled equations of motion for 9 and p. The 
interaction between these two quantities leads to a renor- 
malization of the longitudinal phason velocity. 

The first step in the functional integral formalism is to 
introduce auxiliary fields describing spin and charge fluc- 
tuations. The main technical difficulty is then to recover 
the mean-field (or Hartree-Fock) solution in a saddle- 
point approximation, while maintaining rotational in- 
variance in spin space which is a necessary condition 
for obtaining the spin-wave mode. To overcome this 
difficulty, we introduce a space and time fluctuating 
spin-quantization axis, foilpwing a method introduced by 
SchulzcB and Weng et a/E3 for the isotropic 2D Hubbard 
model. Note that in quasi-lD systems, the distinction 
between right and left moving electrons allows one to 
write the Hamiltonian in a rotationally invariant form 
which is well suited for the calculation of the spin-wave 
mode if we focuses only on the 2kp particle- hole (Peierls) 
channel.ES However, the concomitant consideration of the 
Landau channel (long- wavelength charge fluctuations), 
which is at the heart of our approach, does require the 
introduction of a fluctuating spin-quantization axis. 

The organization of the paper is as follows. The effec- 
tive action of the system is derived in Sec. |l|. We first 
introduce bosonic fields describing charge and spin fluc- 
tuations, and the fluctuating spin-quantization axis. We 
take special care to introduce the physical charge density 
field p. T he standard mean-field theory is recovered in 



calculate the spin stiffness, the spin-wave velocity, and 
the coupling constant of the NLcrM. 

The sliding mode is studied in Sec. IV. We obtain 



Sec. HA, while fluctuations are studied in Sec. [IB. The 



latter are most conveniently computed by performing a 
chiral rotation of the Fermion fields. The corresponding 
Jacobian (the so-called chiral anomaly) is calculated in 
Sec. lilBl . 

The effective action governing the dynamics of the 
spin-wave mode is shown to be a NLcrM together with 
a topological Berry phase term in Sec. HI. We explicitly 



coupled equations of motion for the phase 9 of the SDW 
condensate and the charge fluctuations p. By integrating 
out one of these fields, we obtain the effective action as a 
functional of either 9 or p. We also calculate the conduc- 
tivity and the density-density correlation function from 
the effective action. 

We do not consider long-range Coulomb interaction, 
which would require taking into account normal electrons 
that are thermally excited above the gap. The latter are 
indeed expected-to play a crucial role in the screening of 
the interaction.LJ Note however that the Coulomb inter- 
action affects neither the spin-wave mode nor the trans- 
verse phason mode sampled by optical spectroscopy. 

We consider only the zero-temperature limit, and take 
h — ks ~ i throughout the paper. 



II. EFFECTIVE ACTION 

In the vicinity of the Fermi level, the electron disper- 
sion is well approximated as 

£{kx,ky) = VF{\kx\ - kp) - 2ti,cos{kyb), (2.1) 

where k^ and ky are the electron momenta along and 
across the conducting chains, if, is the transfer integral 
in the transver se d irection, and b the inter-chain dis- 
tance. In Eq. (2.1), the longitudinal dispersion is lin- 
earized around the ID Fermi points ±kp and vp — 



2ata sin(fci?a) is the corresponding Fermi velocity, with 
ta ^ tf, being the transfer integral and a the lattice spac- 
ing along the chains. The linearized dispersion (2.1) sat- 
isfies the property e(k) — — e(k + Q), which corresponds 
to a perfect nesting of the Fermi surface at wave vector 
Q ~ {2kp,TT/b). Actually this property is an artifact 
of the linearization and does not hold for the original 
tight-binding dispersion unless the system is half-filled 
(kp = 7r/2a). Deviations from perfect nesting can be 
taken into account by adding higher harmonics to the 
transverse dispersion —2tb cos{kyb) . For simplicity, we 
shall not consider such terms and restrict ourselves to 
the perfect nesting case. We also assume that the SDW 
is incommensurate with the crystal lattice, so that there 
is no pinning by the lattice. 

Following the standard procedure, we introduce right 
and left Fermionic fields ip+a- and 4'-a- (c is the spin in- 
dex). In terms of the quartet of Fermion fields -0^ = 

('0!|, i/'-ti "04-1 1 "01 1 )? one can write the Hamiltonian of 
the system a.s H — Hq + Hj, where 

Ho = ^ dx \ipi{x)vp{-idxT3 - kp)ipn{x) 

n •' 

-tb X] '4'i{x)i^n+5(x) , 
5=±1 



Ht = U 



dx i^L {x) Ip'l I {x)ll}ni {x)i^n\ (x), (2.2) 



'^ni 



U is the on-site Coulomb interaction strength and n the 
chain index. Here and in the foUowing, r^ (/u = 1,2,3) 
and a^, [v — x, y, z) are 2x2 Pauh matrices acting on 
left/right and spin indices of the Fermionic fields respec- 
tively. The product T^ai, is to be understood as direct 
product of the matrices r^ and ct^ , and any single matrix 
r^ or (Ty as direct product of that matrix with the unit 
matrix. 

Introducing the charge and spin density fields 



Pc = ipH, Ps = V'Vzl/', 

Pc+ ^ llj^T^lp, Ps+ = i'^T^CTzIp, 



(2.3) 



where r± = (Ti±iT2)/2, we rewrite the interaction Hamil- 
tonian as 

Hi ^ j"^ dx [{pcnir, x))"^ - {p,n{T,x)f 



+2(Pc+n(T,2;)Pc+n(T,a;) 

-P*s+ni'r^x)ps+n{T,x))]. 



(2.4) 



In the Matsubara formalism, the partition function 
of the system can be written as a functional inte- 
gral over anticommuting Grassmann variables, Z = 
j Dil)'^Dil)e-^^'^^^'^\ with the action 



S= I dT 
'o 



J2 J dxiji{T,x)drMr,x) + H[^p\^P] 



(2.5) 



where (3 = 1/T is the inverse temperature. The limit 
T ^ is to be taken at the end of the calculations. 

At this point, it is customary to introduce auxiliary 
fields for the spin and charge fluctuations via a Hubbard- 
Stratoppidtch (HS) transformation. However, as noted 
earlieitJ'Ej in the context of systems described by the 
isotropic 2D Hubbard model, such a procedure imme- 
diately leads to loss of spin rotational invariance. The 
reason for this is that in writing down the Hamilto- 
nian, we have made a particular choice (z) for the spin- 
quantization axis of the electrons. Other decompositions 
of Hi in terms of charge and spin fluctuations are pos- 
sible. They are all equivalent as far as the partition 
function is calculat ed e xactly. The reason for choosing 
the decomposition (2.4) is that it allows to recover the 
Hartrce-Fock solution at the saddle poi nt le vel within the 
functional integral formalism (see Sec. II A). 

In order to maintain spin rotational invariance, one 
should consider the spin-quantization axis to be a pri- 
ori arbitrary and integrate over all possible directions in 
the partition function. This is done in practice by intro- 
ducing a new field (/)„ (t, x) which is related to the old 
Fermionic field ipn{T,x) through an unitary SU(2)/U(1) 
rotation matrix Rn{T,x), i.e. ip = Rcj). The rotation 



matrix R satisfies RazR^ = a n, where n is an unit vec- 
tor field which gives the direction of spin-quantization 
axis at point (x, n) and time r. The new field (f) has its 
spin-quantization axis along the local n vec tor. T he i n- 
teraction term in the action, Sj [Eqs. ( ^.4| ) and (2.5)], 
is invariant under this transformation while the unper- 
turbed part, 5*0, becomes 

So = y2 drdx\cbiiT,x)[dr^Co 

+VF{-idxT3 - kp) - VFTzCx\(j)n{T,x) 

-h Y. <t>i{r,x)e-'iT'''''y<'^cl,n+s{r,x)]. (2.6) 



<5=±1 



We have introduced the fields C^ given by 



Co = -R^'drR, 
Cr -= iR^drR, 



Jnb 



Cydy 



RLR 



'71-^-5 • 



(2.7) 



The fields C^ are SU(2) gauge fields and the U(l) gauge 
freedom here corresponds to an arbitrary rotation about 
the z axis, which does not change the state of the sys- 
tem. The partition function now contains an additional 
integral over all C^ or equivalently n field configurations. 
The C^ fields thus contain information about the spin ex- 
citations of the system. In writing Eqs. (2.6) and (pj), 
we have considered only long-wavelength fluctuations of 
the spin-quantization axis, so that the matrix R acts like 
the unit matrix with respect to the left/right indices of 
the Fermionic fields. 

The 2kF charge fluctuations play no role in the SDW 
phase, so that we ignore the term ex p*_|_„pc-(-n in the in- 
teraction Hamiltonian Hj. We thus write the interaction 
term in the action as 

5/ = — V / drdx [{pcn{T,xyf - 2p*_^^{T,x)ps+{T,x)] 

(2.8) 

Note that we have kept explicitly the long-wavelength 
charge, fiuctuations since they couple to the phason 
modeJj In principle, one should also retain the long- 
wavelength spin fluctuations which couple to the spin- 
wave mode. This coupling renormalipes the spin-wave 
velocity by the usual Stoner factor.ld In practice, it 
seems difficult within our formalism to treat the long- 
wavelength spin fluctuations in a way that preserves spin 
rotation invariance. For this reason, we shall ignore them 
in the following. 

We introduce two HS fields, p^^ (real) and A (com- 
plex), corresponding to charge and spin density fluctua- 
tions, respectively. The action of the system can then be 
written a.s S = Sq + Sj, where Sq is given by Eq. (]^.(\), 
and Si by 



r0 , 

^ [^n (t, x)ps+n [t, x) + C.c] 

+ ^[(pf(r,x)f + 2|A„(r,x)|2]}. (2.9) 

Note that the HS field p^^ introduced here is not the 
physical charge-density field, but its conjugate. This can 
be easily checked by varying the action S with respect to 
p . Following Palo et al.^ we now introduce the phys- 
ical charge density field p by decoupling the quadratic 
term in p^ (r, x) by means of a HS transformation. This 
leads to the action 



Si = 



dr dx 



-ipcn{T,x)p^ (t,x) 



U 



IKi'^^ x)Ps+ «(t, x) + c.c] + j(p„(t, x)f 



+ -\^^{t,x)\^ +ip'^'^{T,x)pn{T,x) 



(2.10) 



It can be seen, by varying the action with respect to the 
HS field p^^ , that p is indeed the physical charge density 
field. 

Finally, we introduce a weak external electromagnetic 
field A^ in the action in a gauge-invariant manner. This 
does not change Si while Sq becomes 



So 



x 



V / dTdx\(j>liT, 
-tb Yl <?^Ji('r'2;)e"'J" 

(S=±l 



dr - eAo - Co 



Cx)Tz-kp) (t>n{T,x) 



(Cy+eAy)dy 



(l3n+s{T,x)j. 



(2.11) 

Note that, as an artifact of the linearization of the dis- 
persion relation along the chain, 5*0 does not contain the 
diamagnetic term ex A^. However, as we shall see, this 
term is important for recovering a gauge-invariant effec- 
tive action. We shall recover this term sta rting from the 
original tight-binding Hamiltonian in Sec. HB, 



A. Mean-field theory 



Cfi — Af^ — 0, 
p^^ and A = Aq exp(iQ • r), where Q 



The standard mean-field theory is recovered from a 
saddle-point approximation with n = z 
p = po, p^^ 

is the nesting wave-vector and r = {x,nb). We take Aq 
to be real, without any loss of generality. Varying the 
action with respect to the auxiliary fields p, p and A, 
we obtain the mean-field equations 



Po = --JJPO^ = {(t)li'r,x)(l)niT,x))MF, 



u 



Ao = ^e "^■''(0]j(T,a::)r_tTz0„(r, x))mf- 



(2.12) 



It is clear from the above equations that the mean field po 
is the density of particles in the system and has the same 
value in both the metallic and the SDW ground state. 
As a result, — «Po ^ can be absorbed as a trivial shift in 
the chemical potential of the system [see Eq. ( |2.1C| )]. The 
equation for Aq is the usual mean-field equation for the 
SDW order parameter. 

The mean-field propagator can now be obtained from 
the mean-field Fermionic action 



^MF = 4EE K(-«,k+ Q),0t_^(^.„,k)) g^'M'7:''l'^^ 



^n k,a- 



(/)_cr(Jti^„,k) 



(2.13) 



where we have dropped additive contributions to the 
free energy due to the mean- fields po, p^^ and Ag. Here 
Q~^ is the inverse propagator given by 



g;;' 



iujn - e+{k + Q) sgn(tT)Ao 
sgn(cr)Ao iun-e^ik) 



(2.14) 



where sgn(cr) = -l-(— ) for a =T(i) and e±(k) = vp{±kx — 
kp) — 2ti,cos{kyb) is the dispersion relation for the right 
and left Fermions. 

The mean-field propagator is obtained by inverting 
Ga^- It is given by 

^ _ f G+a{i(^7i,'i^+Q) F+cr(itJ„,k+ Q) \ (0T^\ 

^"^y F_,(ic^„,k) G_,(*c.„,k) )' ^^■^^> 

where 

G±a{iuJn,'k) = -((/)±o.(ia;„,k)(/)^^(iw„,k)) 



e±(k) 



4 



Ag' 



F±c,{iujn,k) = -(0±o-(itj„,k)0|!p^(iw„,k=FQ)) 
sgn(cr)Ao 



'^^ + 4 + Ag' 



(2.16) 



and we have used the relation e_|_(k + Q) = — e_(k) in 
obtaining the above result. 



B. Fluctuations 



We do not consider amplitude fiuctuations, since they 
are gapped and decouple from the sliding and spin- wave 
modes in the long- wavelength limit .111 We therefore write 
the auxiliary fields as 



„HS 



P = Po + Sp, 



(2.17) 



where Sp, Sp^^ , and 9 represent small fluctuations of the 
fields about their mean-field values. The action can then 
be written a,s S = Sq + Sj, with 



^o = E 



dr dx ■ 



[(j>i{T,x)[dr^eAo-C^-i5p'^'' 



+VF {{-idx ~ eA^ - C2:)t3 - kp) 



-Ao e 



,-j(Q-> 



tb ^ 4>i{T,x)e 

fi 

dr dx 



-{n + S)h 



(C: 



<f>n (t, x) 



u 



n 

^(po + Spn{r,x)j 



^{Spn{T,x)f+lSp^^{T,x) 



(2.18) 



The next step to obtain the effective action of the aux- 
iliary fields Sp, Sp^^ and 9 is to integrate out the Fermion 
fields. This is most conveniently done by first introducing 
a new field (j) related to the field 4> through an unitary 
chiral transformation: 



U, 



chiral V 



jT3e/2i' 



(2.19) 



This transformation leaves the interaction part of the ac- 
tion S'j invariant while 5*0 is given by 



^o = E 



dr dx 



{^l^'nHr,^ 



dr 



A tot 

^0 



+VF {-idxTs 



A! 



-Ao(e-'^-'T^a^ -fee 
-tb E <^n^('^'2:)e"*J' 



'/'«(T,a;) 



'"^'''KUy 



<5=±1 



(2.20) 



where the gauge field A*°* is given by 



<' = eAo 



-dx9 + iSp^^ + Co 



Al°' = (eA, - ^—dr9 + Cx)t3, 

2VF 

4°* = eAy - ^dy9T3 + Cy. 



(2.21) 



The chiral transformation therefore eliminates the phas e 
9 of the order parameter A. The action Sq [Eq. ( 2.20 )] 
therefore acquires a simple form. The tradeoff, however, 
is that the Fermions are now subjected to an effective po- 
tential yl*°* which contains the derivatives of the phase 
field 9. Notice that in contrast to the superconducting 



case, the gradients of the phase couple to the external 
electromagnetic field in a different manner. This is a con- 
sequence of the different broken symmetries in the two 
cases. Since we assume that the fluctuations of the order 
parameter (we have taken ^mf = 0) are small, the gra- 
dients of 9 must be small. As a result, one can carry out 
a perturbative expansion of the action in 9^0 or equiva- 
lently in yl*°' . The chiral transformation also produces a 
non-trivial Jacobian J, which yields the additional con- 
tribution Sj = — In J to the action. We shall come back 
to the origin of this non trivial Jacobia n and the method 
of its calculation in more detail in Sec. p^IB 2 . 

We are now in a position to write the action in 
terms of the mean-field action S'mf [Eq. (2.13)] and 
the part involving fluctuations. The action then reads 



S = S'o + S'j + S J, wheve 

Sq ~ Smf 



E 



dr dx 



{^'JiT,x)[Al°'+VFA'f]^'^{T,x) 



-tb E 'Pni'r,x)[(^ 



-^11 



(n + S)b 



Al°'dy 



<5=±1 



0l+5(T-,a;)|. 



(2.22) 



Using the fact that A*°* is a slowly varying weak 

field, we expand the factor exp(— i/^|^ A*°*dy) — 1 
in Sq in powers of A*°*. The terms in the expansion 
are written in a symmetric way with respect to n and 
n + S. Here, we retain terms up to the quadratic order in 
yl*°* . The corresponding contribution to the action reads 






linear 



S^'"", where 



S: 



linear 



~2~ ^ 

n,S=± 



dr dx 



0:t(r,x)(4°*(T,x) 



s 



dia 



+Al°^+s{T,x)jS(b,,+siT,x), 

s 2 , 
+A^yn+8{T,X)^ (J),,+s{t,x). 



E J dTdxc^'Jir,x)[Al°^,ir,x) 



(2.23) 



The linear term Sy'^'^'^'^ corresponds to the coupling of 
A*°* to the paramagnetic current in the transverse di- 
rection, while Sy^'^, which is quadratic in A*°*, is the 
diamagnetic contribution. If we had started with the 
original tight-binding Hamiltonian in the x direction, we 
would have come up, apart from the linear paramagnetic 
term, with a similar diamagnetic term in the x direction 
given by 

,2 rl^ 

ir'- 
,. /o (2-)= 

x\eAxn{Px,T) +Ca;n{Px,T)\'^ 

X [cos(fca;a) -t- cos((fca; - Px)a)](j)'^{kx,T), (2.24) 



*- - 2 ^io "'^ (2^)2 '^"^^^' 



where ta is the hopping parameter along the chain. S^^'^, 
which is not obtained if one uses a hnear dispersion law 
from the very beginning, has to be included in the action 
in order to maintain gauge invariance. 

It is convenient to introduce the charge and spin cur- 
rents for the Fcrmions: 

^ S=±l 



and 5"^'^^ = Sf"" + S^'"" is given by Eqs. (f^) and (|^ 



Integrating out the Fermions, we obtain to quadratic 
the effective action 



order in the fields A^°^ 



C//2 



(2.29) 



MF 



where (• • •)mf means that the average is taken with re- 
spect to the mean- field action ^mf- 

The evaluation of (S"')mf is trivial. Only the term 
involving Jq contributes. We thus obtain 



5 = -^ E ^i^lsr3cl>'n + '/'LV3C5), (2.25) (5")mf = -Po [ drd^ (^5p^^ 



s=±i 

where the index a runs over 0,x,y, z, a is the unit ma- 
trix, jjj is the ^ th component of the charge current, and 
j" for a ^ give different components of the spin cur- 
rent. The current g is almost same as the charge current 
in the transverse direction, except that it is chiral, i.e. it 
has opposite sign for the left and right moving Fermions. 
This current is introduced for notational convenience, as 
we shall see later. 

Furthermore, since the C^ fields are SU(2) gauge fields, 
it is possible to write them in terms of the a matrices, 
namely, C^ = A'^a^, where the index fi runs over 0,x,y 
and the index 1/ over x,y, z. As a result, the fields ^4*°* 
can be expressed as 



/Jtot _ 4O I Al^ 



eAo - ^9a;( 



(2.26) 



where the expr ession s for A'^ can be deduced from 
Eqs. ( I2.21D and (|2.26| ) to be 



A°o = eAo - -j-d^e + iSp 

Al = eA, - -^dr9, 
2vf 



HS 






eAy 



(2.27) 



The action [Eq. ( |2.22|) ] can be conveniently expressed in 
terms of the charge and spin currents as 5 = Smf + S" + 
S'^'^ + S'j + Sj, where 

(2.28) 



(2.30) 



Here, we have taken the continuum limit at the end of 
the calculation and replaced the sum over the chains by 
an in tegra l in th e y d irection. It can be easily seen from 
Eqs. ( |2.18| ) and ( ^.30D that the first term in {S")^^ can- 
cels the term which is linear in Sp in the expression of 
Sj. Thus, we finally obtain 



iS"] 



MF 



-Po 



I dTd\(eAo-^d,9), (2.31) 



and 



S'j^ 



dTffr 



U 



(Spy + iSp^'^sp 



(2.32) 



It may seem unphysical at first sight that the scalar 
potential in the effective action couples to the con- 
stant mean-field density pQ and not to the full density 
P = Pa + Sp, as one would intuitively expect. How- 
ever, one should bear in mind that one stiil has to in- 
tegrate over 6p^^ . As shown by Palo et alSB in the con- 
text of superconducting systems, by redefining the field 
6p^^ -^ Sp^^ - {eAo - ^dxO), we immediately get the 
coefficient of the scalar potential to be the full density p 
and not po. 

In the next three sections, we evaluate the diamagnetic 
term (S"^"')mf, the contribution Sj arising from the chi- 
ral anomaly, and (S"'^)mf- 



1. Diamagnetic contribution 

In this section, we calculate the contribution of t he dia - 
magnetic term to the effective action. From Eq. ( ^.23 ) 
one can easily obtain in Fourier space 



/rtdia\ "-o" ^ 

{by JMF - -J^a 



2 



(3^ 



d'^p 

(2^ 



e'^\Ay{ip„,p) 



x,y,z 



\A;{tPn,p)\' + ^mPn,Pr 



- y. 



■/3 



a— ±,(T,a;„ 






[cos(fcy6) +cos((fcy -py)&)]Gacr(icj„,k), 



(2.33) 



where we use A*(mp,,p) = Ay{—ipn, — p)... for real fields and G±a- is the mean-field propagator given by (2.16). In 
the limit if, ^ AqJ^ we can ignore the effect of the gap and replace G±cr by its value in the metallic phase. To first 
order in Py , we then have 



- y. 



P 



Q — ±,(T,a;„ 



(Pk I P^h"^ 

—--[COS{kyb) + COS{{ky ~ Py)b)]Go.,{lUn,\^) = 4iV ( ) U ^ ^ 



2tt 



cos{kyb) 



X / deQ{—e + 2ti,cos{kyb)), 



2N{0)vj 



, pI^' 



plb"^ 



(2.34) 



where v±^ — ^2t\h'^ is the velocity of the Fermions in the transverse direction, iV(0) = Xj-KVpb is the density of 
states per spin, and O is the step function. Note that only the states near the Fermi level (|e| < 2tb) contr ibute to 
the diama gnetic term (S'^''')mf 
Eq. (|2.33D, we finally get 



so that the use of a linearized dispersion law is justified. Substituting Eq. (2.34) in 



(OMF = A^(0)t'i^E/ (01 (l-^j (e2|^.(*Pn,p)P+ E |A^(*p„,p)P + ^|0(*p„,p)N . (2.35) 



The calculation of the diamagnetic term in the x direction should be done using the original tight-binding dispersion 
law, since it involves electronic states deep in the Fermi sea which forbids the use of a linearized dispersion law. From 
Eq. ( |2.24[ ), we deduce 



/ rtdiaN 

\px /MF 



V-n. \ V 



^\^ 



_fk_ 

(2^)2 



[cos(fc2;a) +cos{{ka: - Px)a)]GaiiUJn,'k), 



(2.36) 



where G^"'^(«a;„, k) = iojn + 2taC0s{kxa) + 2tbC0s{kyb) -\- ji [^ being the Fermi level). Again, we neglect the effect of 
the gap Aq. Neglecting corrections of order tb/ta, we have 



-y 

R ^ 



(fk 
P^J (2^ 

which yields 



[cos(fcj,a) +COii{{k^ -px)a)]Ga[il^n:'^) 



■nb 



''Z" dk 



J /a 



27r 



- cos(fc:ra)6(2ia cos(fc^a) + /i) = -y|iV(0), (2.37) 



{b^ )mf = 



^(o)^^E/ (0^ Ki^-(*p»'P)i' + E i^5;(^p-p 



(2.38) 



r 



We do not consider the term of order g^ , which is con- 
sistent with the linearized dispersion law used in the rest 
of the calculation. 



2. Chiral anomaly 



ter pr oduc es not only a change of the gauge fields A*°* 
[Eqs. ( ^.26 ) and ( 2.27 )] , but also changes the ground state 
of the system. This non-perturbative effect shows up in 
the Jacobian of the chiral transformation. Chiral anoma- 
lies haveJaeea|ki|iasvn for a long time in the context of DW 
systems. ta"E3'c3'Ea Our method of calculation is similar to 
that of Ref. M. 



In this section, we calculate the actio n Sj due to the 
Jacobian of the chiral transformation ( ^.19 ). The lat- 



The chiral transformation changes the local density of 
particle in the ground state (the total particle number re- 



C/chiral[<56'(T,r)] =62 



— „5-r3<5e(T,r) 



(2.39) 



maining unchanged). Since the particle density couples 
to the gauge field Aq, this yields an additional contribu- 
tion, Sj, to the action. 

Let us first calculate the density change 6p{T, r) due to which changes the phase of the order parameter from 
an infinitesimal chiral transformation, 6{t, r) to d{T, r) — S9{t, r): 
I 

6piT, r) = lim (0t(r + 6T,r + Sr) (ul,^,,JSe{T + 5r, r + <5r)]C/ehirai[<^^(r, r)] - l) 0(t, r)\ , (2.40) 



where "»e,use the point spitting regularization scheme. A regularization is necessary to properly calculate the particle 
density.cJ To lowest order in 69, we obtain 



Sp{T,r) = -- lim [5e{T + ST,r + Sr)-S9{T,r)]{(j)^T + ST,r + Sr)T3(t>{T,r)) 

2 (5r,(5T— ^0 

= -^ lim [de{T + dT,r + 5r)-Se{T,r)]y"aGaa{-6T,-6r), 

2 Sr.Sr^O ^ — ^ 



(2.41) 



where Gaa is the mean-field propagator defined in Sec. IIA.Lj Introducing the Fourier transforms S6{ipn,p) and 
GQo-(icj„, k), Sp is rewritten as 



6p{T, r) 
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d^p 



<fk 



(27r)2 J (27r) 



50{ipn, p)e*P-^-*^"^ ^ a[G„,(ic^„ + Jp„, k + p) - G„,(^c^„, k)], (2.42) 



after a trivial shift of integration variables, and in the limit 6t = 5y = 0. Performing the sum over cj^, we obtain 

5p{T, r) = -^J2j^J^ ^^('P-^ p)e^P'--^"^ J2 "[""(k + p) - n„(k)], (2.43) 



where ^^(k) = 0(fci? — ak^ + {2ti, / v f) c,os{kyb)) is the particle occupation number. [Note that nQ.(k) has the same 
value in the metallic and SDW phases at T — Q] Using 



9/ \ / *?/ 

Ua^k- + p) — "■Q(k) = —iapx -\ Py sin(fcj,6) \5[kp — ak^ H cos{kyb) 

vp / V vp 



(2.44) 



r 



we finally obtain, to lowest order in p, 
i -^-^ f (Pp 






p 



pje{ipn,p)e'P 



r — ipnT 



Sj = -vpN{0) / d'rdr d^e{eAo + i<5p"^) 



(2.47) 



1 



= —rdx 

TTO 



(T,r). 



(2.45) where we have used l/nb = vpN{0). 

In deriving the action Sj, we have implicitly assumed 
By integrating out the Fermions, we have shown that that the chiral transformation only produces a change 



the ground-st ate de nsity of Fermions couples to the gauge 
field Aq [Eq. ( 2.30| )]. Since an infinitesimal chiral trans- 



in the density of the grou nd state. Following analogous 
steps leading to Eq. (2.4E), it can be shown that the chi- 



formation changes the density by S^it also produces the ^.^i transformation does not change the spin density and 



contribution to the effective actionEJ 



SSj = - 



d'^rdrdpA^ 



-— / d'^rdrdjefeAo + i6p^^ - —dj 
TTbJo V 2 . 



(2.46) 



Taking a variational integral of Eq. ( p.46| ) , we obtain the 
chiral anomaly contribution to the effective action 



the spin or charge currents, so that the change in the 
density is the only non-perturbative effect of the chiral 
transformation. 



3. {S"''h 



In this section, we calculate {S"^)mf- Using ( 2.2g ), we 
deduce 



{S )mf 



~ -^ ^ ' a/3 Ml' 



/3 



P71 

0^x,y,z 0.x,y 



Yl Y ^"iiPn,P)^3^g{iPn,P)iPy6{-ipn,-p) + ^\0{ipn,P)fngg{iPn,P) 



(2.48) 



where n^.„^.^(ip„,p), Ujo,g{ipn,p) and ngg(ip„,p) are 
the current-current correlators defined by 

^i°jf («Pn' P) = Om (*P»' P)J^{-iPn, -P))mF ' 
%g(*Pn>P) = (i;^(«Pn,p)g(-«Pn,-p))j^p, 

ngg(ip„,p) = {9{iPn,p)9{-iPn, -P)>MF ' (2-49) 

and the currents j'!^{ipn, p) and g{ipn, p) are the Fourier 
transfo rms o f the charge, spin and chiral currents defined 
inEq. (|^). 

In the long- wavelength limit, we need to compute the 
effective action to second order in a gradient expansion. 



For symmetry reasons, the correlators IIx,;, o,z, 11 



n,x,« , and Hi^n vanish, and IIix^x 






We also 



to first order and II^o^o to second order. As shown in 



the Appendix, we thus obtain 



find that none of the spin currents j'f^'^'^ couple to the 
charge currents j^ or g, which simply states that the slid- 
ing and spin-wave modes decouple in the long-wavelength 
limit. Since the SU(2) gauge fields are first order in 
gradient {A" oc da-n"), we need H-, „' to zeroth order. 

This also holds for the correlator Ilgg , since it couples to 
[dyO)'^ . Based on similar reasoning, we must obtain Ilgjo 



n^oAiPn,p)=N{0) 



vIpI + v\pI 



3AI 



mojo{ip„,p) = 2N{0)vl 1 



6Ag 



2 I 1 i^y" 



n,o,o(zp„,p) = 2A^(0K 1- 



n.o,.o(zp„,p) = A^(o)4-^"^^ 



6A2 



n-oo(ip„,p) = N{0)vl 



9 '^PnPy 



%io" (iP", P) = Ujyjy {ipn, p) = 2iV(0), 

%j-(*P„,p)-2A^(0)4, 

^&ii^Pn:P)^2N{0)vl. 



(2.50) 



Separating charge and spin contributions, we have 

-(5"2)mf/2 = ^2tpi„ + ^2^^hason, whcre 



5*; 



iff 

2 spin 



N{O)^J2J-0j2[\^o(.Wn,p)\' + \Ayi^p„,p 



Ceff 

'^2 phasoii 



-^(0)^Y|^{\'M^Pn,p) 



)|2 + 4|A^(*p„, p)|2 + «^|A^(jp„, p)|^ 



HS/- ^\\2 i '^FPx+^l_Py \ 1\A I- „m2 2 



+ i(5p"=(ip„,p)|^ 



6Ag 



-^ e^\ASPu.p)\vp 1 



6A§ 



+e^|A,(zp„,p)|M(l-^ 



P« 



6A2 / ^ 4l^(*P"'P)I^Pn - d{ipn,p)evFPnA.j:{-ipn-,-p) 



lp„ 



3A2 



+iX^ [e^o(iPn,p) +i(5p"^(ip„,p)] [u|.p:re^:r(-ip„, -p) + w5^pj^eylj^(-ip„, -p)] 

I 



This completes our derivation of the effective action. 
In the next sections, we shall use this effective action to 
obtain the spin-wave and the phason modes. 



A. NLcrM 



(2.51) 



The contribution to the spin part S'^pj^ of the effective 
action comes from the ter ms i nvolvi ng th e gauge fields 
A^^v,^^ From Eqs. ( ^^ , ( |I3^ ) and (|]5|), we find 



III. SPIN- WAVE MODE 



^eff 



-vl\Al{ip^,p)f - vl\A;,{ip^,p)W 



(3.1) 



The c ontri bution from the diamagnetic terms [Eqs. ( 2.35| ) 
and ( 2.38 )] cancels the terms proportional to A^ in 



Eq. (2.51), so that we are finally left with an effec- 
tive action which depends only on Af^ and A^. This 
cancelation is expected since our effective action has 
to be gauge invariant. As mentioned earlier, for the 
SU(2) fields C^ a gauge transformation corresponds to 
an arbitrary rotation about the z axis, and does not 
change the state of the system. Such a rotation chants 
A^ -^ A^^ + 1/2 9^A, where A is the rotation angle.tJ 
Thus to be gauge-invariant, the effective action cannot 
depend on the A^ field. The diamagnetic terms in the 
effective action are absolutely crucial for this cancelation 
and hence for obtaining a gauge- invariant effective action 
for the spin-wave mode. 
Using the identitiesEj 



tr[Cl) - tria.Co)' = Af + Af = -i^, 



(3.2) 



we express the effective action [Eq. (3.1)] in terms of gra- 
dients of the n field: 



for |pn I , "Wf IPx I , i'± by I ^ Ao, the NLcrM describes the 
behavior of the model only i n th e low-energy long- 
wavelength limit. Thus, Eq. ( |3.3| ) should be supple- 
mented with cut-offs Aq ^ Aq in energy space and 
Ax ^ Aq/vp and Aj, ^ Ao/wj_ in momentum space. 

The actioD (3.3) can be expressed in a more conven- 
tional formLj as 

5.fi„ = 11 d'rdr [x(a.n)2 ^ ^^^(5^„)2 + p^^i^dyuf] , 

(3.4) 

where we have introduced the parameters x = A^(0)/2, 
p^x = 7V(0)w|/2 and p^y = N{0)v']_/2. Note that 
these values are only approximate. In order to ob- 
tain the correct low-energy long- wavelength behavior, 
it would be necessary to integrate out all high-energy 
short-wavelength fluctuations. The latter are expected 
to renormalize the bare values of %, p^^ and psy quoted 
above. 

To see the physical interpretation of the parameter 
X, we note that x = ^f^j^iWu = 0,p = 0)/4 (i/ = 

x,v)- The correlators II,-^,-^ and II.B.y can be linked 

. " " .... •'o.-?o 

to the transverse spin susceptibilities in the mean-field 

state, i.e. H, 



Jo -^0 



= 4(S';j,(y)S'j;(y))MF, where S^ 



s: 



pi = ^l' d'-dr [(9.n)^ -t- 4(5.n)^ + .i(9,n)^] . 

(3.3) 

The effective action for the spins therefore turns out 
to be an anisotropic NLaM, as expected. Since S^^^^ 
has been obtained within a gradient expansion, valid 



(f)^a^4)/2. We see that x = {SxSx)mf = {SySy)MF 
and hence corresponds to the transverse static uniform 
spin-susceptibility of the system in the mean- field ground 
state. 

The spin stiffness coefficients psx and psy , on the other 



hand , come from the diamagnetic terms [Eqs. ( 2.35 ) and 
( p. 38 )], which are themselves related to the average ki- 
netic energy. To see this point more clearly, let us con- 
sider the diamagnetic term in the y direction. We have 
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^(0KE^E/7^K(^^"'P)I'-^E^E/7^I^^(^^«'^ E 7|^2cos(M)G..(^^„,k) 



%LJ„ .a.a 






i2n) 



KiiPn^pW 



2 1 y 



(3.5) 



where {Ky)MF is the average kinetic energy in the y direction per unit area in the mean- field state. But we can also 
write the diamagnetic term as 



N{0)v 



x-.y 
,^r^ 



^E/i?.Kc...p)i-fr-v(i 



(3.6) 



h\Ky) 



MF 



Comparing Eqs. (^) and ( p.6[ ), we immediately see 

that the spin stiffness parameter psy is related to the av- ^^ 4 

erase kinetic energy in the y direction. A similar result ^m i • . • <■ , i • i 

, , , . , P f f, 1, , ihe velocities of the spin- wave mode, Vx 



can be obtained for psx, and we finally get 

a^{Kx)MF 



and Vy = [psy/xY^ are then given by 



(3.7) 



Psx 



-a'iKx) 



MF 



4x 



Vp, 
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-b"^ {Ky) MF 

4x 



1/2 



(3.8) 



We find that the longitudinal velocity equals the Fermi 
velocity, because we have neglected the coupling of the 
spin-wave mode with the long-wavelength fluctuations. 
Taking into account the latter would lead to the renor- 
mahzed spin-wave velocity v^ — Wi?(l — [/iV(0))^/^.El 

Finally, we note that we can carry out an appropriate 
rescaling of lengths to obtain an isotropic NLcrM. Rescal- 
ing the lengths as x' = xjvx, y' = v/vy and r' = r, we 
obtain 



"feff 
-'spin 



\fPsxPsy 



d^r'dr' J2 (^^n) 



(3.9) 



The cut-off A' ~ Aq (in reciprocal space) is now isotropic. 
The di mens ionless coupling constant of the isotropic 
NLcrM (3^) equals Ao/ipsxPsy)^^^- It is proportional 
to Aq oc e"^/^''^'-"-' and heiice approaches zero in the 
weak-coupling limit {U ^ 0)Hj 



B. Berry phase term 

In general, in addition to the usual dynamical terms, 
the action describing antiferromagnetic spin systems in 
one or two dimensions is expected to have a topological 
Berry phase term.E3 Such a term was derived for antifer- 
romagnetic systems descijibed by the 2D isotropic Hub- 
bard model by Wen et alt^ We now derive this term for 
our action. 

The Berry phase term comes from the term propor- 
tional to Aq in (5"')mf- However, to get a non-zero re- 
sult, it is not sufhcient to retain only the long- wavelength 
part of the R matrix, since the Berry phase term precisely 
results from the 2fci? oscillating part. Instead of consid- 
ering the more general structure of the R matrix, we go 
back to the original tight-binding dispersion law, which 
is a more direct way for obtaining the Berry phase term. 
The term proportional to Aq in (5"')mf then reads 



T-)az4>)MF in the continuous form ulatio n, we obtain, us- 
ing the mean- field equations [Eq. ( 2.12| )] , 



^1 



Berry 



d^r dT Al{(j)'^ [to 
4Ao 



;mf 



U 



■ / d^rdr AqCOb{Q ■ r). 



(3.11) 



In general, it is not possible to express the field Aq in 
terms of the n vector, because the former is a gauge- 
dependent quantity. Nevertheless, it is still possible ta 
express the variation of the Aq field in terms of n asEj 
Aq — —i/2JQ dr]n{ri)-{dnn{T])xdTn{i])), where 77 is an ex- 
ternal parameter varying continuously from to 1, n(l) 
is the physical n field and n(0) = 0. The action can then 
be written as 



2A f 
•S'sorry = «-^ / d'^r dr cos(Q • r) 



X d77n(77)-[a„n(77)x9,n(r;)]. (3.12) 
Jo 

Another equivalent way of writing the Berry term in the 
effective action is to express the gauge field Aq in terms 
of the polar and azimuthal angles a and (j) of the n vec- 
tor. Using Aq = —i/2 (1— cos(a(r, r))dr(f>{T, r), the Berpy. 
term in the action then assumes the more familiar formEEl 



^Bcriy 



/ d^rdr cos(Q • r) 



.2Ao 
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U 

X [1 — cos(a(r, r))]dr(j){T, r). 



(3.13) 



For 2D and 3D antiferromagnets, it is generally be- 
lieved tliat the presence of the Berry phase term plays 
no role.ES The dynamics of the spin-wave mode is there- 
fore determined by the NLctM [Eq. (3.4)] derived in the 
previous section. 



IV. PHASON MODE 



■Secrry = - / d^r dr A^ ((?!)V 



z(P/MF, 



(3.10) The contribution to the phason eff ective action com es 
from the diamagnctic terms [Eqs. ( 2.35 ) and ( 2.38 )], 
where {4>''(Tz4')mf is the spin density in the mean-field the chiral term Sj [Eq. (b^47)] , the interaction term 5*/ 
state. Noting that (<^V2(/))mf becomes {4)'^{tq + tj^ + [Eq. ( ^.32 )] and 5'2phason [Eq. ( |2.5l|) ]: 
I 



O^x.y 



1 f d^p ( 1 

'S'phason = -3^ j T^ I XI 2 [(e^M(«P»'P) +«'^P"^(«Pn'P)^M0)^P>'(*Pn,P)(eAi.(-ip„,-p) +?:(5p"^(-ip„,-p)5i.o)] 
0,x,y 

+ X! [(e^A'C^Pn'P) + iSp^^ {ipn,p)Sf,o)Ffj,{ipn,p)9{-ipn, -p)] + -\9 {ipn , p)\'^ D {ip„ , p) 

(4.1) 



+iSp{ipn,p)Sp^^{~ipn,-p) + -j\Sp{ipn,p)\'^j, 
where the polarization tensor P^^, and the coupling coefficients F^ and D, are given by 
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3A^ Pn 



(4.2) 



The con tributio n fro m the diamagnetic terms in 
Eqs. (2.3f;) and (2.38) exactly cancels the gauge- 
noninvariant terms in S'2phason- In fact, it can be ex- 
plicitly checked that the polarization tensor P^^ is trans- 
verse. Also, we note that the phase field 9 couples to 
the (gauge- invariant) electric field Ex = —Pn-^x +iPx-Ao- 
This, together with the transverse polarization tensor 
Pf_i^, ensures that if we integrate out the fields Sp, Sp^^ 
and 9 to obtain an effective action iS"'^[Ap], it will be 
gauge-invariant. Since the coupling term between 9 and 
Aq comes from Sj, it turns out that the contribution 
from the Jacobian of the chiral anomaly is crucial for 
obtaining a gauge-invariant effective action. 



A. Equations of motion 

We consider the electromagnetic field as being the ep. 
ternal field, i.e. we neglect the phason-polariton mode.Lj 
Varying the effective action with respect to the fields 9, 
p and p^^ , we obtain the three coupled equations 

D{ip„,p)9{ipn,p) + ^ Ff,{ipn,p){eAi,{ip„,p) 



vfN{0)px9{ip„, p) + iSp{ip„, p) = 0. (4.6) 

Thus, we obtain the well-known relation between charge 
and phase fiuctuations: 

SpiT, r) - VFNiO)dx9{T, r) = \dx9{T, r). (4.7) 

TTO 

From (^^), we finally obtain the (decoupled) equations 
of motion 

..2 q2 2 c,2\ 



-|-i(5p"==(ip„,p)(5^o) =0, 



(4.3) 



Foi-ipn, -p)9{ipn, p) + Sp{ipn, p) 



0,x,y 



^ Pot,{ipn,p){eA^,{ip,i,p) +i5p^^{ip„,p)Sf,o) = 0, 



—Sp{ipn,p)+iSp^^{ipn,p) = 0. 



(4.4) 
(4.5) 



Using Eq. (4.f;), we can eliminate the HS field Sp^^. 
This gives two coupled equations for the physical fields 
6p and 9. Subs tituting the expressions for D, P^^ and F^ 
from Eq. (4.2), and retaining terms up to second order 
in frequency and momenta, we get 

{pl+vlpl+vlpl)9{ipr„p) 

-ivFPxUdp{ipn,p) + 2evFEx{ip„,p) = 0, 



{dt - vi,X - vid;)9{t, r) = -2evFEx{t, r). 



(df - vl^d^ - vldl)6p{t,r) = -^d,Ex{t,r), (4., 



where we have Wick-rotated back to real time t, and 
■^ph = V F {l+U N {0)y /^ is the renormalized phason veloc- 
ity. This renormalization of the longitudinal phason ve- 
locity has been obtained previously from a more conven- 
tional approach based on Green function calculations .El 



B. Effective action S[e,Af,] 

In this section, we derive the effective action S'^^[9, A^j] 
for the phase mode. 

If we try to integrate out the HS fields 6p and 6p^^ in 

[Sp, 6p^^ , 9, Afj], we immediately face the problem 

of inverting the coefficient Pqo of the quadratic term in 



Coff 
phason 



Sp^^, since Pg^^ is singular at Pn,P = 0. Nevertheless, 
we note from the results of the last section, that this 
term does not contribute to the equations of motion up 
to quadratic order in external momenta and frequency. 
We can therefore ignore this term while integrating out 
the field Sp^^ . As we shall see, the effective action so ob- 
tained, although not exact, reproduces the correct equa- 
tions of motion for the phase field 9. 

Integration of the 6p^^ fields within the above-stated 
approximation then gives a factor of J^ p^i^P ~ 
ipx9 /lib) and we are left with the partition function 



Z ^ / DSpD9DAf,S{Sp -Px9)e 



±nn]^-S''''lSp.A,,0]_ 



nb 



(4.9) 



with the action 
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0,x,y 



1 p jz ^T^^y 1 

O.x.y ^ ^ 

+ Y, [eA^{ipn,p)F^,{ipn,p)e{-ipn,-p)] + -\e{ipn,p)\^D{ipn,p) + j\Sp{ipn, p)\^ j . (4.10) 



Integrating out 5p, one then gets the effective action 



^'='^[0, A^] = ^11 / 7^ { II 2 M/.(*Pn,P)^M-(«P",p)e^.(-ip„,-p) 

Pn M^ 

+eiV(0)?;i.^^(ip„,p)6l(-ip„,-p) + -|6'(ip„,p)|2i:)'(ip„,p)|, 



(4.11) 



r 



where L''(ip„,p) = N{0){pl + w^ p2 + vlpl)/2 is the 



C. Effective action S[5p,A^] 



effective inverse phason propagator. This clearly demon- 
strates the modification of the phason velocity due to 
the interaction between the phase field 9 and the density 

fluctuation 6 p. It is easy to check, by varying Eq. (4.11) To obtain the effective action S[Sp,A^] for the den- 
with respect to 9, that this effective action reproduces sity fluctuations Sp , we start from the effective action 
the same equation of motion for 9 as Eq. ( [l.6| ) derived in S'^^[Sp,9,Afj] [Eq. ( |4.10| )], and integrate out the phase 
the earlier section. fleld 9. One then obtains the effective action 
I 

S[Sp,Af,] = -gY j T^ { XI 2 Mp(*Pn,p)-PM'^(JPn,p)eAj.(-ip„,-p)] 

-e(Ao(ip„,p) - -^A^(ip„,p)](5p(-ip„,-p) + -\5p{ipn:P)?Xpp\. (4-12) 



where we have introduced the density-density correla- 
tion function 

Xpp(«Pn,p) = (bp{ipn,p)^p{~ipn-,-p)) 



Mt,r)^~—dt9{t,r) 



(4.15) 



2N{0) 



2 2 
VpPx 



pI + ^IhPl + '"IpI ' 



(4.13) 



It is easy to see from the expression of the density-density 
correlation function that the phason does not induce any 
transverse density fluctuations. 

The currents in the longitudinal and transverse direc- 
tions are given by 



jx{iPn,P) 
Jv{-iPn,P) 



ss 



SAx{~ipn,-p) 

ss 



SAy{-ipn,-p) 



A„=0 



A„=0 



epn 
iPx 



Note that the expressions of the currents as a function 
of the phase 9 [Eqs. (4.7) and (4.15)] have been obtained 
to lowest order in a gradient expansion. They have been 
recently generalized to higher order in Ref. EH. 

The current-current correlators can also be ob taine d 
from the effective action. Using Eqs. ( 4.13| ) and ( 4.14 ), 
we find that the current-current correlator Xji^j^ — {jp.ji^) 
can be expressed in terms of the density-density correla- 
tion function as 

e2p2 

Xj.jAiPn,P) = 5^Xpp(JPn,p) 



dp{ipn,p), 



(4.14) 



pI 



= -2N{0)e' 



2 9 



Pn + ^IhPl + '"IpI ' 



X]y]y{iPn,P) = X]^]y{iPn,P) = 0. 



(4.16) 



These expressions agree with the (real time) continuity 
equation dtp + V ■] = 0. Also, we see that there is no cur- 
rent due to the phason mode in the y direction. The con- 
tribution to the current across the_chain comes entirely 
from the quasi-particle excitations.E3 From Eqs. (4.7) and 
(4.14), we obtain (in real time) 



From the expression of the current-current correlator, we 
then obtain the ac conductivity axxii^) as 



0-x£c(w) = Xj^j, (iPn, P = 0) 

Pn 



2N{Q)e^vl ( iV 



-f7rJ(w)), (4.17) 
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where V denotes the principal part. The real part of the 
conductivity therefore satisfies the f-sum rule 



Re[crj;3;(w)]da; = — , 



(4.18) 



where Wp — i^e^vp jhY^"^ is the plasma frequency. The 
phason mode completely exhausts the spectral weight 
and we obtain the well-known resultcl that there is no 
contribution to the longitudinal conductivity from the 
quasi-particle excitations in a clean SDW system. 



V. CONCLUSION 



term. By a suitable length rescaling, one can obtain an 
isotropic NLcrM with an effective dimensionless coupling 
constant Ao/(ps2^Psj,)^' ^, where p^x and psy are the spin- 
stiffness coefficients. The coupling constant, proportional 
to the mean-field gap Aq oc e^'^l^^'S^) ^ is small in the 
weak coupling limit UN{0) ^ 1. 

The sliding mode is governed by an effective La- 
grangian C{9,p) which is a function of two independent 
fields: the phase 9 of the SDW condensate and the charge 
density field p. From the coupled equations satisfied 
by these two variables, we obtain the phason dynamics 
and its contribution to the current-current and density- 
density correlation functions. 



We have derived the effective action for the low-energy 
collective modes of quasi-lD SDW systems. The in- 
troduction of a fluctuating spin-quantization axis in the 
functional integral allows to consider the phason mode, 
the spin- wave mode, and the long-wavelength charge fluc- 
tuations on equal footing. 

We find that the spin-wave mode is governed by an 
isotropic NLcrM together with a topological Berry phase 
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APPENDIX A: 



In this section, we sketch the calculation of the current-current correlators. These correlators are given by 

n,-.,.' (ipn, P) = n^jv (ip„, p) = {j{ipn, P)j'i-ipn, -P))mF ' (^1) 

The currents j^ (or j in the simphfied notation) can be defined using Eq. ( p.2q ) as 



j-(jp„,p) = ^=^ 



d'k +- ^' 



VP^J (27r)2 



y^ y^ '^io-' (^^", ^)Vaaa' {ky, ky + Py)(j)aa{ii^n +ip„,k + p). 



(A2) 



Q aa' 



where we have set the area of the system to unity. The function s Vga a' (ky , ky + Py) are the current operators which 
can be easily read off from the deflnition of the currents in Eq. ( 2.25| ). For example, for the charge current in the x 
direction j^, Vaaa'iky, ky + Py) — vpaS^^i. From (A2), we then deduce 



+- n 



T^jj'{iPn,p) = --^X! / 72^ SIX! [^o.<^<^'iK - Pv,ky)v^^,^{ky,ky - py)Gaa'{i(^n,'k)Gaa{iuJn -ip„,k-p) 

i^n ^ ' a era' 

+Vaaa'iky - Py , ky)v^^, ^{ky +TT/b,ky -py + 7r/b)Faa'iiuJn,'k-)Facr{iuJn - ip„,k- p - ttQ) , (A3) 

where a = =F for a = ± and G±a and F±a- are the mean-field propagators. We evaluate the frequency sum and the 
kx integral using eo,(k — aQ) — — ea(k) and expanding the energy dispersion e^ in powers of external momenta: 



ea(k) - ea(k - p) = e^; -|- €y, 

ey = 2tbbpy sm{kyb) - Ub'^pl cos{kyb) + 0{p^y), 
ex = avFPx- 



(A4) 



Using these relations and Eq. ( 2.16| ) one can evaluate the frequency sum and the kx integral up to second order in p„ 
and ex = Ca: + ey (0(p^, e|,PneT)) to obtain 
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R Z^ I -^ Faa'[l^n,k)Faa[iUJn - iPn,ii- p - aQ) = Sgn(cr)sgn(cr ) ^— [^ UA^ 



13 y^. TT^ '^aa' {iUJn, k)Gaa{i^n " iPn, k - p) = - 



2 \2 

N{0) (1 pl-e\ ^ ipneT 



UAl 



6Ag 



(A5) 



Substituting Eq. (A5) in the expression for the correlator [Eq. (A3)] we get 



njj'(jp„,p) 



A^(0) f dky 
"2^ 



+- Ti 

EE 



'^acrcr'yr^y Pyi '^y)'^ aa' aV^y ■> ^y Py) 



1 pI 



12A§ 



iPn£T 

6A2 



+Vaaa' [ky - Py , ky)v^„, ^{ky + ir/b, ky~py + 7r/6)sgn(cr)sgn(cr') f - 2 + 



1 , Pn + Sq 



12A2 



(A6) 



The dependence on ky comes from ey and Vatrtr' ■ The latter is either ky independent or depends on ky through 
the function / = sin((/cy — Py)b) + sin{kyb) = 2 sin(A:j^6) — Pybcoa{kyb) — p^b"^ sm{kyb)/2 + 0{py). To evaluate the ky 
integral in (|A6|), we therefore need only a limited set of integrals, which, up to order Py, are given by 



(/)fey - {<^v)ky ~ {f^y)ky - (/ iy)ky - 0, 

(/'>fe„-2-^, {eDuy^vlpl, {fey)ky^2hbpy, {feD^y ^ ivlp^ 



(A7) 



where we hav e us ed the sho rt-hand notation (• • ■)k for J dky /{2Tr) ■ ■ ■. 

From Eqs. (A6) and ( A7), it is then a matter of straightforward algebra to obtain the different correlators. As an 
illustrative example we sketch the calculation of IIjOjO . Since the current operator in this case is given by Vacru' (ky, ky — 
Py) = tbb6aa'f{ky,py), we obtain 



njO,o(ip„,p) 



^^£i:*-(/^fi-''" ■'"■•" 



2N{0)vl (l - ^^ 



pI 



6Ag 



6Ag 



6Agr 



(A8) 



where we have used Eq. (A7) to evaluate the ky integrals to obtain the last line and retained only terms up to second 
order in external frequency and momenta. 

The other correlators can be obtained following an exactly similar procedure and in this way, we finally obtain 
Eqs. (2.50) for the charge and spin current-current correlators. 
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